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SOME  CONSTITUTIVE  RESTRICTIONS  IN  PLASTICITY 


P.  M.  Naghdi 

University  of  California,  Berkeley 


ABSTRACT 

An  account  is  given  here  of  some  recent  results  in  the  mechanical  theory  of 
finitely  deformed  elastic-plastic  materials,  which  are  chiefly  concerned  with 
some  constitutive  restrictions.  These  restrictions  are  discussed  both  in  the 
context  of  the  general  theory  and  for  special  materials.  The  development 
presented  includes  also  a discussion  of  an  alternative  formulation  of  plasticity 
theory  relative  to  loading  surfaces  in  strain  space  and  elaborates  on  its  sig- 
nificance, as  well  as  several  of  its  features. 


1.  INTRODUCTION 

Starting  with  the  general  constitutive  equations  of  the  purely  mechanical 
theory  of  elastic-plastic  materials  contained  in  the  work  of  Green  and  Naghdi 
[1,2]  , the  object  of  this  paper  is  to  present  an  account  of  some  aspects  of 
recent  results  in  plasticity  given  by  Naghdi  and  Trapp  [3-6].  In  the  discus- 
sion that  follows,  we  freely  quote  from  these  works  and  for  details  refer  the 
reader  to  the  original  papers  [1-6]. 

After  some  background  material  and  preliminaries  in  Section  2,  the  main 
features  and  significance  of  an  alternative  formulation  of  plasticity  theory 
relative  to  loading  surfaces  in  strain  space  are  described  in  Section  3»  In 
Section  4,  we  first  introduce  a physically  plausible  assumption  concerning  the 
nonnegativeness  of  the  work  done  by  external  forces  on  the  body  in  a closed 
cycle  of  spatially  homogeneous  deformation.  From  this  assumption  and  with  the 
use  of  certain  results  from  Section  3>  we  then  derive  a work  inequality  for 
finitely  deformed  elastic -plastic  materials  which  holds  in  any  smooth,  spatially 

^The  theory  of  elastic-plastic  materials  in  [1,2]  includes  thermal  effects 
and  is  developed  within  the  framework  of  thermo-mechanical  theory.  Although  we 
confine  attention  to  the  purely  mechanical  aspects  of  the  subject,  the  various 
results  given  here  can  be  easily  interpreted  in  the  context  of  the  isothermal 
theory.  The  purely  mechanical  theory  of  elastic-plastic  materials  with  large 
deformations, which  forms  the  starting  point  of  this  paper, corresponds  to  the 
general  developments  in  [1]  prior  to  thermodynamic  restrictions.  In  Section  2, 
we  summarize  the  form  of  the  theory  of  elastic-plastic  materials  given  in 
Section  4 of  [2] . 


homogeneous,  cloved  cycle  in  strain  space.  Next,  with  the  use  of  the  latter 
inequality,  two  local  inequalities  are  obtained  which  place  restrictions  on  the 
constitutive  equation  for  the  plastic  strain  rate.  The  rest  of  the  paper 
(fection  5)  deals  with  constitutive  restrictions  for  special  elastic-plastic 
materials,  including  certain  geometrical  interpretations  and  restrictions  on 
the  stress  response  for  finitely  deformed  ductile  metals.  The  cases  of  rigid- 
plastic  materials  and  elastic-plastic  materials  with  small  deformations  are 
also  discussed  in  Section  5> 

While  the  alternative  formulation  of  plasticity  relative  to  loading 
surfaces  in  strain  space  (Section  3)  is  also  useful  in  the  derivation  of  some 
cf  the  restrictions  discussed  in  Lections  4 and  5,  the  potential  utility  of 
this  formulation  in  numerical  computations  should  not  be  overlooked.  In  con- 
trast to  the  formulation  of  plasticity  relative  to  loading  surfaces  in  stress 
space,  the  constitutive  equations  for  the  plastic  strain  rate  and  the  associated 
loading  criteria  in  the  alternative  formulation  hold  without  change  in  the  case 
of  elastic-pcrfectly  plastic  materials  and  in  the  entire  physically  work- 
hardening range  of  metallic  materials.  Related  additional  remarks  are  made  in 
Lection  3* 


2.  PRELIMINARIES.  GENERAL  BACK®OUIJD. 

Since  many  concepts  concerning  the  mechanical  behavior  of  materials  are 
extensions  of  observations  made  in  simple  tension  or  simple  shear,  it  is  help- 
ful to  recall  briefly  the  familiar  idealized  response  of  a work-hardening 
elastic-plastic  material  in  a one-dimensional  test  in  which  the  deformation  is 
homogeneous.  In  the  absence  of  hysteresis,  the  nature  of  the  response  of  a 
typical  ductile  metallic  material  in  a simple  tension  test  is  depicted  by  the 
curve  OABH  in  the  plane  of  the  conventional  engineering  stress  (rr)  versus 
engineering  strain  (e  of  Fig.  1.  From  the  point  0 to  the  proportional  limit  A 
the  material  Is  linearly  elastic  and,  since  the  deformation  in  this  range  is 
reversible,  unloading  takes  place  along  AO.  For  loading  above  A (regarded  to 
be  coincident  with  the  initial  yield  point)  the  deformation  is  irreversible  and 
the  material  work-hardens  along  ABH.  Unloading  from  a point  B is  assumed  to 
take  place  elastically  along  BD;  the  compressive  yield  limit  at  D is  now  lower 
than  that  at  A (due  to  Bauschinger  effect)  and  BD,  in  general,  is  taken  to  be 
parallel  to  OA.  Reloading  from  a point  such  as  C (or  any  point  on  BD)  proceeds 
along  CB  (since  hysteresis  is  assumed  to  be  absent)  to  the  point  B,  whereupon 
further  loading  proceeds  to  subsequent  yielding  along  the  path  BH. 

The  schematic  diagram  represented  by  the  solid  line  curve  in  Fig.  1 does 
not  prominently  display  the  extended  range  of  the  rising  portion  of  the  n-e 
curve;  however,  it  does  also  exhibit  the  falling  portion  of  the  curve,  as  well 
as  the  point  H (corresponding  to  the  ultimate  strength)  at  which  the  slope  of 
the  stress-strain  curve  vanishes.  The  corresponding  stress-strain  curve  for 
an  elastic-perfectly  plastic  material  (the  dashed  line  OAA'A")  is  also  shown 
in  Fig.  1.  Ordinarily  in  the  development  of  the  theory  of  elastic-plastic 
materials,  regardless  of  whether  the  deformation  is  finite  or  infinitesimal, 
attention  is  confined  to  the  rising  portion  of  the  n-e  curve  (OAB)  along  which 
the  slope  of  the  stress-strain  is  positive  and  the  falling  portion  of  the  curve 
(beyond  H)  is  ignored.  However,  it  is  important  to  note  that  a theory  which  is 
constructed  for  the  work-hardening  range  of  the  metallic  materials  may  not  be 
valid  in  the  vicinity  of  the  ultimate  strength  (or  ultimate  load)  or  in  the 
region  about  which  necking  begins?  We  return  to  this  aspect  of  the  subject  in 
the  next  section. 

Let  the  motion  of  a body  be  referred  to  a fixed  system  of  rectangular 
Cartesian  axes  and  let  the  position  of  a typical  particle  in  the  present  con- 
figuration at  tine  t be  designated  by  x^  = Xj^(X^,t),  where  XA  is  a reference 
position  of  the  particle.  It  is  convenient  to  introduce  the  symmetric  strain 
eKL  by 


The  term  necking  is  used  here  synonymously  with  the  maximum  point  on  the 
rr-e  curve  (see  Fig.  3*20,  p.  87  of  Richards  [7])* 


Fig.  1 An  idealized  diagram  for  a work-hardening  material 
(curve  OABH)  in  the  n-e  plane  exhibiting  a region  of  zero 
slope  and  the  corresponding  stress-strain  curve  for  an 
elastic-perfectly  plastic  material  (the  dashed  lines 
OAA'A").  Also  shown  is  unloading  elastically  from  the 
point  B along  BC. 


where  cKL  are  the  components  of  the  Cauchy-Green  measure  of  deformation  defined 
in  terms  of  the  deformation  gradient  F.j  K relative  to  the  reference  position  and 
5KL  denotes  the  Kronecker  symbol.  Throughout  the  paper,  all  subscripts  take 
the  values  1,2,3  and  the  usual  summation  convention  is  employed  over  repeated 
subscripts . 

The  deformation  function  Xp  in  (1)^  may  be  specified  in  the  form  Xp  = 
xK5Ki+ui’  where  up  = Up^XM,t)  are  relative  finite  displacement  components. 
For  later  reference,  consider  now  a simple  tension  test  from  which  the  standard 
engineering  stress-strain  curve  such  as  the  n-e  curve  of  Fig.  1 can  be  obtained 
Let  the  initial  length  of  the  tensile  specimen  be  £0,  its  final  length  i and 
choose  the  Xp-axis  to  coincide  with'the  longitudinal  axis  of  the  specimen. 
Further,  let  Xp  = Xp  + Up(Xp)  be  the  deformation  function  in  a homogeneous 
extensional  deformation,  where  the  relative  longitudinal  displacement  Up(Xp) = 

[ (i-20)/i0]Xp,  and  let  e denote  the  resulting  Lagrangian  strain  [corresponding 
to  the  component  epp  in  (l)p].  Then,  with  the  use  of  (l)p,  the  engineering 
strain  e defined  by 

du.  l-l  du  . 

e=^=-T^  > F,  = 1 + — t = = 1 + e (2) 


We  further  recall  the  relationship  between  the  non-symmetric  Piola- 
Kirchhoff  stress  rr. K,  the  symmetric  Piola-Kirchhoff  stress  sKT  and  the  symmetric 


Cauchy  stress  j , namely 

niK  " ?iLSLK  ’ SLK  SKL  * °ij  = [det(FkM)J  F'iKFjLEKL  * ^ 

The  stresses  niK  and  o- , are  sometimes  referred  to  in  the  literature  as  the 
engineering  (or  nominal;  stress  and  the  true  stress , respectively.  Let  tt  and  s 
refer,  respectively,  to  the  engineering  stress  [corresponding  to  tt,.  in  (4),] 
and  to  the  component  s-q  of  the  symmetric  Piola-Kirchhoff  stress.  Then,  with 
reference  to  the  one-dimensional  hcmogeneous  extensional  deformation  referred 
to  above,  by  (4)^  and  (2)2  we  have  the  relationship 

rr  = (l+c)s  . (5) 


It  is  easily  seen  from  (3)  and  (5)  that  if  the  deformation  is  infinitesimal  so 
that  e«l,  then  (3)  reduces  to  e = e to  the  order  of  approximation  considered 
and  the  distinction  among  tt,s  and  the  corresponding  true  stress  (i.e.,  the 
component  in  (4),)  -disappears. 

We  summarize  now  the  principal  results  for  the  mechanical  theory  of 
finitely  deformed  elastic-plastic  continua  contained  in  the  more  general  thermo- 
dynamical work  of  Green  and  Naghdi  [1,2].  Thus,  in  addition  to  the  strain  e^ 
defined  by  (1)^>  at  each  point  of  the  continuum  we  admit  the  existence  of 
(i)  a plastic  strain3  specified  by  a symmetric  second  order  tensor  eBr;  (ii)  a 
measure  of  work-hardening  specified  by  a scalar  function  K;  and  a scalar-valued 
function  f -- called  the  yield  or  loading  function  — which  depends  on  the  vari- 
ables s^j,ePj,  and  K and  which  is  such  that  the  equation 


f(sMN’eM’K)  = 0 


(6) 


for  fixed  values  of  ej^j  and  K represents  a hypersurface  in  the  six-dimensional 
Euclidean  space  of  the  symmetric  stress  sKL.  We  assume  the  loading  function  f 
to  be  continuously  differentiable  with  respect  to  its  arguments  and  that  the 
yield  condition  (6)  be  an  orientable  surface  of  dimension  five  in  the  six- 
dimensional stress  space. 

Let  the  strain  eKL  depend  on  the  set  of  variables 


and  write 


U = (sMN’e£N’K) 


(7) 


eKL  eKL^  * 

Then,  assuming  that  the  expression  (8)  is  invertible  for  fixed  values  of  e!L 
and  ic,  we  have 


s 


KL 


(U) 


(9) 


where  for  convenience  we  have  introduced  the  abbreviation 


^ (eMN,eMN,IC)  ‘ (10) 

The  response  of  elastic -plastic  materials,  which  is  developed  relative  to  the 
loading  functions  in  stress  space  and  which  is  confined  to  the  work-hardening 
range  of  the  material,  is  characterized  by  the  constitutive  assumption  (9),  as 
well  as  those  for  the  rate  of  work-hardening  and  the  plastic  strain  rate, 
namely 


and 


K hKLeKL 


(n) 


^Although  we  assume  here  that  e^L  is  symmetric,  the  developments  that  fol- 
low can  be  modified  to  include  the  effect  of  a non-symmetric  plastic  strain 
tensor  included  in  [2]. 


i 

1 


4. 


when  f < 0 


when  f = 0 and  f < 0 


when  f = 0 and  f = 0 


f when  f = 0 and  f > 0 


where  a superposed  dot  stands  for  material  time  derivative  with  respect  to  t 
holding  XK  fixed,  and  h]^  are  functions  of  the  set  of  variables  (7)> 

f is  defined  by 


and  where  the  partial  derivative  Bf/dSj^j  Stands  for  the  symmetric  form 
2(31/92^1  + af/dsj.jjj) . The  four  conditions  involving  f and  r on  the  right-hand 
sides  of  the  four  equations  in  (12)  are  called  the  loading  criteria.  Using  the 
conventional  terminology,  these  four  conditions  in  the  order  listed  correspond 
to  an  elastic  state,  unloaaing  from  an  elastic-plastic  state,  neutral  loading 
and  loading  from  an  elastic-plastic  state. 

The  constitutive  equations  (12)  are  not,  in  general,  valid  for  the  special 
case  of  an  elastic-perfectly  plastic  material  or  in  regions  where  the  material 
behavior  corresponds  to  that  at  ultimate  strength.  In  fact,  in  the  context  of 
the  formulation  of  plasticity  relative  to  loading  surfaces  in  stress  space,  the 
theory  of  elastic-perfectly  plastic  materials  is  usually  developed  separately. 
As  noted  by  Green  and  Naghdi  [1],  the  constitutive  equations  for  an  elastic- 
perfectly  plastic  material  in  the  presence  of  finite  deformation  may  be  deduced 
as  a limiting  case  of  the  above  results  after  allowance  is  made  also  for  a 
modification  of  the  loading  criteria.  Thus,  if  k is  always  constant  ( = K0  say) 
and  the  yield  function  f is  independent  of  ej^,  the  loading  surface  is  always 
stationary  and  may  be  represented  as 


Moreover,  instead  of  (12),  for  an  elastic-perfectly  plastic  material  we  have 


when  f,  < 0 


when  f. 


when  f. 


where 


pPQpFQ.  ° MN  ™ 

While  the  criteria  for  an  elastic  state  and  for  unloading  in  (15 )q  2 31,6  Eijnilar 
to  those  in  (12 )j_  o > it  should  be  remembered  that  neutral  loading  Ao  longer 
exists  in  the  elastic-perfectly  plastic  case  and  that  the  criterion  for  loading 
in  (15)3  does  not  involve  a condition  similar  to  f >0  in  (12)^. 

The  form  of  the  nonlinear  theory  of  elastic-plastic  materials  outlined 
above  represents  the  material  description  of  the  theory.  It  is  developed  rela- 
tive to  loading  surfaces  in  the  six-dimensional  stress  space  and  is  intended  to 
be  valid  in  the  work-hardening  range  of  the  material  corresponding  to  the  range 
in  which  the  stress  is  monotonically  increasing  function  of  strain  in  one- 
dimension.  It  is  perhaps  worth  recalling  that  the  development  of  the  constitu- 
tive results  which,  in  particular,  leads  to  (11 ) and  (12)  begins  with  the 
assumption  that  the  response  functions  for  e^L  and  K are  linear  in  the  stress 
rate  sKL  and  that  both  eKL  and  k are  independent  of  the  particular  time  scale 
used  to  compute  the  rate  of  change.  Further,  both  eKL  and  eB^  are  regarded  as 
independent  kinematic  variables  and  only  the  tensor  eKL  is  related  to  the 
deformation  gradient: 


Although  the  stress  tensor  which  naturally  occurs  in  the  equations  of 
motion  of  the  material  description  is  the  non-symmetric  Piola-Kirchhoff  stress, 
the  variables  which  occur  in  the  constitutive  equations  are  the  variables  b 
which  include  the  symmetric  Piola-Kirchhoff  stress.  For  example,  instead  of 
(6),  we  may  begin  by  considering  a yield  function  which  depends  on^  tt^k.  But 
this  stress  is  related  by  (4)^  to  the  symmetric  stress  tensor  sK^  and  the 
deformation  gradient,  so  that  the  yield  function  nay  be  expressed  as  another 
function  of  sK^  and  F^.  Moreover,  the  yield  function  must  remain  unaltered 
when  the  continuum  is  subjected  to  superposed  rigid  body  motions  and  it  follows 
that  the  yield  function  can  be  expressed  in  terms  of  a different  function  of 
Skl  and  eKL.  In  view  of  the  constitutive  assumption  (8)  for  the  strain  tensor, 
the  yield  function  can  be  reduced  further  to  the  form  in  (6). 


3.  AN  ALTERNATIVE  FORMULATION  OF  PLASTICITY  IN  STRAIN  SPACE 


In  order  to  motivate  the  significance  of  an  alternative  formulation  of 
plasticity  theory  relative  to  loading  surfaces  in  strain  space,  we  now  compare 
the  response  of  a typical  ductile  metal  in  the  one-dimensional  homogeneous 
extensional  deformation  considered  previously  with  the  corresponding  prediction 
of  the  theory  of  elastic-plastic  materials  summarized  in  Section  2.  As  noted 
with  reference  to  Fig.  1,  the  uniaxial  stress-strain  curve  resulting  from  a 
simple  tension  test  is  a plot  of  engineering  stress  n as  a function  of  engi- 
neering strain  e , while  the  corresponding  response  curve  according  to  the 
theory  of  Section  2 is  a plot  of  the  symmetric  Piola-Kirchhoff  stress  s as  a 
function  of  the  Lagrangian  strain  e.  Keeping  this  in  mind  and  recalling  the 
relations  (3)  and  (5),  it  can  then  be  shown  that  [3] 


CURVE  A 


CURVE  B 


STRAIN 


STRAIN,  in  / in.  x 10 


Fig.  2 Conventional  engineering  stress-engineering  strain  diagram 
(the  n-e  curve  in  the  lower  left-hand  corner)  for  7075  aluminum  alloy 
together  with  the  curve  A representing  an  enlarged  portion  of  the  rt-e 
curve  and  the  curve  B representing  the  corresponding  portion  of  the 
s-e  curve.  The  approximate  coordinates  of  points  M and  M'  are 
(4.6  xlO-2,  84xlo3)  and  (4.7  xlO-2,  8o.3xlo3),  respectively. 


yThe  yield  function  may  also  be  assumed  to  depend  on  F 
is  not  essential  to  the  argument  that  follows. 


The  n-e  stress-strain  curve  shown  in  the  lower  left-hand  side  of  Fig.  2 
corresponds  to  the  experimental  data  for  7075  alluminum  alloy  given  in  Richards 
C 7 ; p.  137 » Fig.  3»55].  Curv"  A in  Fig.  2 is  an  enlargement  of  the  indicated 
portion  of  the  n-e  curve  and  curve  B represents  the  corresponding  portion  of 
the  s-e  curve  calculated  with  the  use  of  (3),  (5)  and  (17).  Several  features 
of  the  curves  A and  B should  be  noted  and  compared: 

(a)  The  slope  drr/de  in  (17 )_  is  always  positive  since  n/(l+e)>0  and 
hence  (17)2  represents  the  slope  on  the  rising  portion  of  the  n-e  curve  prior 
to  necking  or  ultimate  strength; 

(bj The  s-e  curve  will  always  have  a zero  slope  at  a point  within  the 
work-hardening  range  of  the  n-e  curve  and  the  maximum  value  of°  s occurs  at  N' 
directly  below  N on  curve  A and  to  the  left  of  the  point  M at  which  n is  a 
maximum ; 

(c)  The  slope  ds/de<0  fior  the  range  of  values  of  the  stress  between  N' 
and  M',  which  corresponds  to  f < 0 and  e^=0  in  (12),  and  thus  is  contrary  to 
the  physical  observation  that  the  increment  of  plastic  strain^  is  not  zero 
between  N and  M ^or  between  N'  and  M'); 

(d)  Since  f=0  at  N',  as  this  point  is  approached  in  (12 )l  must 

become  unbounded  and  this  renders  the  value  of  the  plastic  strain  rate  inde- 
terminate at  N'. 

In  contrast  to  the  singular  behavior  at  N'  on  curve  B,  the  corresponding 
point  N on  curve  A is  in  the  physically  work-hardening  range  of  the  stress- 
strain  curve  and  also  in  the  region  prior  to  necking  or  initiation  of  insta- 
bility. Moreover,  (12)4  will  be  nearly  singular  over  a range  of  stress  around 
N'  and  the  loading  criteria  in  (12)3  4,  which  involve  the  stress  rate,  fail  to 
hold  at  all  points  beyond  N'.  This  means  that  the  plastic  strain  rate  is 
extremely  sensitive  to  the  stress  rate  and  hence  a representation  of  the  form 
(12)4  may  not  be  reliable  in  any  explicit  numerical  calculations  in  the  vicinity 
of  a point  such  as  N'.  It  is  worth  noting  also  that  the  particular  features 
associated  with  the  zero  slope  of  the  stress-strain  curve  in  a simple  tension 
test  do  not  arise  in  a one-dimensional  homogeneous  compressive  deformation. 

This  is  because  in  the  compression  test  the  cross-sectional  area  of  the  speci- 
men is  continuously  increasing  (instead  of  decreasing  in  the  case  of  a simple 
tension  test)  and  hence  there  are  no  points  of  zero  slope;  see,  in  this  con- 
nection, Richards  [J;  p.  150,  Fig.  4.3} • 

Motivated  by  the  above  background  and  in  order  to  obtain  an  appropriate 
expression  for  the  plastic  strain  rate  which  would  be  valid  in  the  full  range 
of  plastic  deformation,  we  summarize  briefly  here  the  main  aspects  of  a dif- 
ferent formulation  of  plasticity  in  the  form  given  by  Naghdi  and  Trapp  [3l- 
The  loading  criteria  in  this  formulation,  as  well  as  all  constitutive  equations 
corresponding  to  (ll)  and  (12),  are  expressed  as  a linear  function  of  e^,  in 
addition  to  their  dependence  on  the  variables  (10).  Thus,  we  admit  the  exis- 
tence of  a loading  function  g (e^jej^jjic)  such  that  the  equation 

g(eMN’eMN’K)  = ° (18) 

for  fixed  values  of  ej^  and  K represents  a hypersurface  of  dimension  five  in 
the  six-dimensional  Euclidean  space  of  the  symmetric  strain  e^.  Alternatively, 
(18)  can  be  obtained  from  (6)  after  substituting  (9)  for  the  argument  s^  in 
the  function  f.  Again  as  in  the  theory  of  Section  2,  we  assume  that  the  rate 
of  work-hardening  and  the  plastic  strain  rate  are  independent  of  the  particular 
time  scale  used  to  calculate  the  rate  of  change.  Then,  with  the  assumption 
that  the  plastic  strain  rate  depends  linearly  on  e^  and  employing  a procedure 
similar  to  that  used  in  [1],  the  constitutive  equations  for  k and  ej^  can  be 
expressed  in  the  forms 


The  maximum  value  of  s on  curve  B is  calculated  by  means  of  (3),  (5)  and 
(17)  from  the  same  experimental  data  but  plotted  in  the  s-e  plane. 

7 

Recall  that  since  the  deformation  is  homogeneous  the  increment  of  plastic 
strain  corresponds  to  the  rate  of  plastic  strain. 


— 


I 


7. 


mKLeKL 


^PKLg 


when  g < 0 , 

A 

when  g = 0 and  g < 0 , 

A 

when  g = 0 and  g = 0 , 

A 

when  g = 0 and  g > 0 


where  and  m^  are  functions  of  the  set  of  variables  (10)  and  g stands  for 

-4^-  (21) 

Supplementary  to  the  results  (20),  it  is  also  shown  in  [3]  that  when  the  stress 
loading  criteria  imply  one  of  the  conditions  in  (12)]_  2 a (corresponding  to  an 
elastic  state,  unloading  and  neutral  loading,  respectively),  the  strain  loading 
criteria  in  (20)-^  g 3 inply  precisely  the  same  conditions  relative  to  strain 
space.  As  far  as ’loading  is  concerned,  no  general  conclusion  can  be  reached 
regarding  the  correspondence  or  equivalence  of  g>0  and  f >0,  but  we  emphasize 
that  the  constitutive  equation  (20)^  holds  in  the  full  range  of  plastic  deforma- 
tion. Moreover,  in  contrast  to  the  modified  form  of  the  loading  criteria  for 
the  elastic-perfectly  plastic  case  in  (15 ) as  compared  with  those  in  (12),  the 
loading  criteria  in  (20)  of  the  alternative  formulation  hold  without  change  for 
elastic-perfectly  plastic  materials. 

A discussion  of  comparison  of  the  initial  and  subsequent  yield  surfaces  in 
strain  space  with  corresponding  results  in  stress  spaces  is  illustrated  in  [3] 
with  the  use  of  two  well-known  yield  functions,  namely  those  of  von  Mises  and 
Tresca.  In  particular,  it  is  shown  that  the  initial  yield  surfaces  of  von  Mises 
and  Tresca  have  exactly  the  same  geometrical  relationship  to  each  other  in 
strain  space  as  they  have  in  stress  space.  For  further  discussion,  we  refer 
the  reader  to  [3,  Sec.  5]« 


4.  RESTRICTIONS  ON  CONSTITUTIVE  EQUATIONS  FOR  PLASTIC  STRAIN  RATE. 

We  recall  here  certain  restrictions  derived  by  Naghdi  and  Trapp  [4]  from 
consideration  of  nonnegative  external  work  in  a cycle  of  homogeneous  deforma- 
tion and  then  elaborate  on  the  nature  of  the  results. 

Consider  a closed  cycle  of  spatially  homogeneous  motion  in  the  closed  time 
interval  [t-^^],  (t^<t2)V  The  cycle  is  said  to  be  smooth  if  the  time  deriva- 
tives of  displacement,  strain  and  associated  kinematic  quantities  are  continuous 
in  C t_i  , t2 ] and  assume  the  same  values  for  each  material  point  at  times  t^  and 
t^.  We  designate  such  a smooth  spatially  homogeneous  closed  cycle  of  deforma- 
tion by  C(t2?t2)  and  recall  the  following  work  assumption:  The  external  work 
done  on  the  body  by  surface  tractions  and  by  body  forces  in  any  smooth  spa- 
tially  homogeneous  closed  cycle  is  nonnegative,  i.e., 

PlL  "ikVi"*!  >»Vidvlit  > “ <22> 

t-,  ft 


for  all  cycles  C(t-L,t2).  In  (22),  R0  :.s  the  region  of  space  occupied  by  the 
body  in  its  reference  configuration,  dRQ  is  the  closed  boundary  surface  of  R0, 
p0  is  the  mass  density  in  the  reference  configuration,  v^  are  the  components  of 


Recall  that  a homogeneous  motion  is  one  whose  deformation  gradient  is 
independent  of  the  material  coordinates  so  that,  in  a spatially  homogeneous 
motion,  the  strain  tensor  e^p,  is  a function  of  time  only,  for  a closed  spa- 
tially homogeneous  cycle  in  the  closed  time  interval  ft^^],  the  displacement 
x^  and  the  strain  e^  assume  the  same  values  at  times  t^  and  t^. 
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the  velocity,  are  the  components  of  the  body  force  per  unit  mass,  NK  are  the 
components  of  the  outward  unit  normal  to  dftQ,  represent  the  components  of 

the  stress  vector  acting  on  the  body  in  the  present  configuration  and  measured 
per  unit  area  in  the  reference  configuration,  and  dA  and  dV  refer  to  elements 
of  area  and  volume  in  the  reference  configuration. 

We  now  proceed  to  derive  a certain  inequality,  which  involves  the  stress 
power,  from  the  assumption  (22)  and  the  equations  of  motion 


iK,K  Ko“i  ‘ Vi  ’ 

where  a comma  denotes  partial  differentiation  with  respect  to  Xr.  Consider  the 
inner  product  of  (23)  with  v^,  integrate  the  resulting  expression  over  the 
region  R0  and  use  the  divergence  theorem  to  obtain 


where  in  writing  the  integrand  of  the  last  integral  we  have  used  (4) 
and  the  identity 


iK  i,K  it  LK  i , K KL  KL  ' ' ' 

Next,  integrate  (24)  in  the  time  interval  [t^jtg]  and  combine  the  result  with 
the  assumption  (22)  to  arrive  at  the  inequality 

[2{£\  5 P^Vj v . dV  + I*  sne  dV}dt  a 0 . (2< 


Since  the  closed  cycle  of  deformation  ^(t-ptg)  is  smooth,  the  velocity  (and 
therefore  the  kinetic  energy)  has  the  same  value  at  t = tj_  and  t = t2>  so  that 
the  first  term  in  (26)  vanishes  for  every  cycle  cCt^jtg).  Further,  since  Sj^ 
and  e^L  are  functions  of  time  only  and  the  volume  J dV>0,  (26)  or  equiva- 
lently the  assumption  (22)  can  be  reduced  to9  140 


The  inequality  (27)  derived  from  (22),  is  valid  for  any  smooth  closed  homogene' 
ous  strain  cycle.  Now,  let  ejj^  be  any  strain  inside  the  loading  surface  in 
strain  space  so  that 

g(eMN’eMrK)  < ° ’ (28 

and  consider  a smooth  strain  cycle  in  strain  space  associated  with  C(t-j_,t2) 
beginning  and  ending  at  eSr.  For  this  cycle,  (27)  can  be  written  as 


er  integrating  the  last  result  by  parts  and  remembering  that  e5Qj(t^_) 
(tg ) = e®^,  we  obtain 


The  inequality  (29)  holds  for  any  smooth  closed  strain  cycle  beginning  and 
ending  at  ej^  inside  the  loading  surface  (18)  in  the  strain  space.  It  may  be 
contrasted  with  a similar  inequality  over  a closed  stress  cycle  in  the  stress 
space,  namely 


^An  inequality  similar  in  form  to  (27)  but  in  the  context  of  small  deforma 
tion  is  the  starting  point  of  Il’iushin's  discussion  in  [9]  and  is  referred  to 
by  him  as  the  'postulate  of  plasticity. ' 


(30) 


(6KL-SKL)eKLdt50 


where  s^..  is  any  stress  inside  the  loading  surface  (6).  It  should  be  noted 
that  in  the  context  of  finite  deformation,  (30)  cannot  be  deduced  from  the 
assumption  (22).  However,  the  linearized  version  of  (30)  for  small  deforma- 
tions corresponds  to  Drucker's  Postulate  [10];  in  this  connection  see  also  [11]. 

With  the  use  of  (29)  and  by  considering  a special  sequence  of  spatially 
homogeneous  closed  strain  cycles,  two  inequalities  are  derived  in  [4]  which 
hold  during  loading.  One  of  these  inequalities  can  be  used  to  prove  that  a 
tensor  gKL  defined  by  % A 

g =!!kl-p  (3D 

gKL  p MN  3*  * 1 

°em 

is  directed  along  the  normal  to  the  loading  surface,  i.e., 

a =-v  d£i  . yiQ 


gKL  “ * Y 


where  y is  a nonnegative  scalar  function  of  the  variables  (10),  and  the  second 
inequality  is  given  by 

gKL4KL  5 0 * (33) 

Now  for  certain  purposes,  it  is  more  convenient  to  express  the  stress  response 
in  terms  of  a different  but  equivalent  set  of  kinematic  variables  (e^-ejk, 
e?T,t().  Thus,  for  the  stress  constitutive  equations,  we  may  also  write  (see 

[67) 

6KL  = “tt(eMN-eMN>*MN»,C)  (34) 

A 10 

and  also  assume  that  the  response  function  s in  (9)  are  such  that 

a A kb 


With  the  use  of  (34),  (35)  and  the  chain  rule  for  differentiation,  it  is 
shown  in  [4]  that  the  results  (32)  and  (33)  can  be  referred  to  stress  space 
and  expressed  in  the  forms 


' ®KLMNe 
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r%  pP  p SO 
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where  the  fourth  order  tensor  functions  are  defined  by 

_ !!kl  + !!a  h i 

HOMN  3s  „ L n 3k  MNj 
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Thi^  symmetry  restriction  is  tantamount  to  assuming  that  the  response 
functions  s^,  are  derivable  from  a potential.  In  fact,  if  the  constitutive 
equations  (6)  to  (12)  are  viewed  in  the  context  of  the  isothermal  theory  of 
elastic-plastic  materials  in  which  s^L  are  expressed  in  terms  of  the  partial 
derivatives  of  e^,  then  (35)  is  automatically  fulfilled. 


It  should  be  noted  that  the  restrictions  (36 -_>  hold  in  all  notions,  even 
though  they  have  been  derived  iron  consideration  of  homogeneous  strain  cycles. 
Indeed,  since  the  response  functions  which  occur  in  (36)2  2 depend  only  on  the 
local  kinematic  variables  (10)  and  not  their  spatial  gradients,  it  is  sufficient 
to  consider  spatially  homogeneous  motions  in  order  to  obtain  constitutive 
results  which  would  be  valid  in  all  motions. 

A _The  coefficients  -Bp and  in  (37)  involve  fhe  response  functions_ 

®KL*SFCL  ant*  ®KL*  once  s-r  is  known,  its  inverse  e^j  and  the  functions  Cj^ 

are  also  determined.  Hence,  according  to  (36)0  the  plastic  strain  rate  can  be 
determined  from  the  knowledge  of  th<=  loading  function  f,  the  work-hardening 
response  function  hj^  and  the  stress  response.  It  is  clear  from  (36)1  that 
only  in  special  cases,  e.g.,  when  are  assumed  to  be  indepedent  of  ej^j  and 
K , * ne  plastic  strain  rate  will  be  directed  along  the  normal  to  the  yield 
surface  in  stress  space.  A further  discussion  on  this  will  be  postponed  until 
Section  5- 


5.  SOME  SPECIAL  RESULTS  AND  SPECIAL  ELASTIC-PLASTIC  MATERIAL. 

We  consider  in  this  section  some  special  elastic-plastic  materials  and 
special  cases  of  the  results  discussed  in  Section  4,  including  a geometrical 
interpretation  of  the  inequality  (29)  in  one-dimension  and  restrictions  on  the 
stress  response  of  ductile  metals.  Among  the  various  cases  considered  are 
several  special  materials  for  which  stronger  restrictions  than  those  indicated 
in  (36)2  2 can  be  found. 


(A)  A Geometrical  Interpretation  of  the  Inequality  (29) 
and  its  Consequences 

In  order  to  gain  sane  insight  into  the  nature  of  the  inequality  (29) 5 which  is 
a cons  quence  of  the  assumption  (22),  in  this  subsection  we  consider  a geo- 
metrical interpretation  of  a one-dimensional  version  of  (29)  and  compare  this 
with  a corresponding  inequality  over  a stress  cycle.  To  this  end,  again  let  s 
and  e stand  for  the  one-dimensional  components  of  the  stress  s^  and  the  strain 
e^  in  a homogeneous  finite  extensional  deformation  of  a ductile  material  and 
let  (s°,e°)  refer  to  an  elastic  state  prior  to  any  additional  plastic  deforma- 
tion. Then,  in  one-dimension  and  over  a strain  cycle,  (29)  can  be  written 
as 

t 

I ^(e-e°)sdt  S 0 or  | (e-e°)ds  5 0 , (38) 

Jt1  JC(e) 

where  C(e)  in  (38)g  refers  to  a closed  strain  cycle  such  as  the  cycle  ABCDF  in 
Fig.  3 and  where  an  increment  of  stress  ds  (since  the  deformation  is  homogeneous) 
corresponds  to  sdt.  The  inequalities  in  (38)  are  obtained  from  (29).  They 
should  be  contrasted  with  similar  inequalities  in  one-dimension  and  over  a 
stress  cycle  which  follow  from  (30),  i.e. 

["  ^(s-s°)edt  ? 0 or  f (s-s°)de  5 0 , (39) 

\ JC(s) 

where  C(s)  in  (39)2  refers  to  a closed  stress  cycle  such  as  the  cycle  ABCD  in 
Fig.  3 and  where  an  increment  of  strain  de  (since  the  deformation  is  homogeneous) 
corresponds  to  edt. 

With  reference  to  the  s-e  diagram  in  Fig.  3s  we  observe  that  on  the  rising 
portion  of  the  one-dimensional  stress-strain  curve  the  negative  of  the  left- 
hand  side  of  (38)2  represents  the  area  enclosed  by  ABCDFA  (cross-hatched 
vertically  in  Fig.  3)  and  the  left-hand  side  of  (39 )o  represents  the  area 
enclosed  by  ABCDA  (cross-hatched  horizontally  in  Fig.  3)-  Parallel  statements 

^The  line-rized  version  of  (39)  corresponds  to  a statement  of  Drucker's 
postulate  [10]  in  one-dimension.  A discussion  of  some  aspects  of  cons equen^.s 
of  (38)  and  (39)  > but  limited  to  small  deformation,  is  contained  also  in  a 
paper  by  Palmer,  Maier  and  Drucker  [12]. 
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concerning  the  interpretations  of  the  left-hand  sides  of  (38)2  and  (39)?  hold 
also  for  the  falling  portion  of  the  s-e  curve  over  the  strain  cycle  a'BT'D'F 
and  over  the  stress  cycle  A'B'C'D'  in  Fig.  3- 


Fig.  3 An  idealized  one-dimensional  mechanical  response  of  a work- 
hardening material  (curve  OABCB'C'M)  in  a homogeneous  extensional 
deformation  plotted  in  the  s-e  plane.  The  symbols  (s°,e°),  (s'°,e'°) 
refer  to  an  elastic  state  or  prior  to  any  additional  plastic  deforma- 
tion. Also  shown  is  unloading  elastically  from  various  points  on 
both  the  rising  and  the  falling  portions  of  the  curve,  as  well  as 
the  closed  strain  cycles  ABCDF  and  A'B'C'D'F'  and  the  corresponding 
closed  stress  cycles  ABCD  and  A'B'C'D'. 


While  both  (38)  and  (39)  can  be  assigned  the  geometrical  interpretations 
throughout  the  entire  hardening  range  of  the  stress-strain  diagram,  the  conse 
quences  of  the  inequality  (39)  cannot  be  valid  on  the  falling  portion  of  the 
s-e  curve.  To  see  this,  consider  the  inequality  {36)?  when  the  coefficients 
CkLMN  are  constants  and  independent  of  the  variables  flO) . Then,  for  an 
isotropic  material  and  in  one-dimension,  it  can  be  easily  shown  that  (36)2 
yields 


where  C is  a material  constant.  At  a typical  point  on  the  fairing  portion  of 
the  s-e  curve,  say  at  the  point  P in  Fig.  3»  we  have 


If  we  also  place  the  usual  restriction  on  the  coefficient  C (corresponding  to 
Young's  modulus  of  elasticity),  we  see  that  (40)  is  satisfied  and  the  work 
assumption  (22)  is  consistent  with  the  physical  behavior  of  ductile  metals  in 
the  falling  portion  of  the  one-dimensional  stress-strain  curve.  In  contrast 
to  the  observation  just  made,  the  inequality  corresponding  to  (40)  resulting 
from  (39)  over  a stress  cycle  is  of  the  form  seP  ? 0 and  this  is  clearly 
violated  over  the  falling  portion  of  the  s-e  curve  where  s<0  and  eP>0  by 


(B)  Restrictions  on  the  Stress  Response  of  a Class  of 
Elastic-Plastic  Materials 

As  indicated  in  [6],  a properly  invariant  representation  can  be  obtained  for 
the  stress  and  the  free  energy  response  functions  in  terms  of  easily  inter- 
pretable kinematic  measures.  Such  representations  can  then  be  used  to  further 
restrict  the  form  of  the  stress  response  function.  With  reference  to  Fig.  1, 
recall  now  that  in  the  case  of  ductile  metals  unloading  from  an  elastic- 
plastic  state  takes  place  linearly  along  a straight  line  such  as  BC,  even  at 
large  (or  moderately  large)  strains.  Making  an  assumption  on  the  basis  of 
this  observation,  an  expression  for  the  stress  response  is  obtained  in  [6] 
which  is  linear  in  certain  kinematic  measure  (essentially  the  elastic  strain 
tensor)  but  with  coefficients  which  are  arbitrary  functions  of  the  plastic 
strain  and  the  work-hardening  parameter.  For  details,  the  reader  is  referred 
to  [6],  where  a discussion  of  a number  of  other  special  cases  can  be  found. 


(C)  Materials  Whose  Stress  Response  is  Linear  in  (e._-e?„T) 

MH  MM 

As  noted  by  Kaghdi  and  Trapp  [5],  in  the  construction  of  the  proofs  in  [3] 
which  led  to  the  results  (36),  an  explicit  use  was  not  made  of  the  stress 
response  in  the  form  (9)  and  only  an  expression  for  the  rate  of  stress  s^  was 
utilized.  Consider  now  a special  class  of  elastic-plastic  materials  whose 
stress  response  may  be  specified  through  the  rate  of  stress  in  the  form 
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where  the  tensor 
conditions  (42)2' 
vanish  and 


KLMN 
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(42) 


coefficients  Lf™  depend  on  ejU  and  K and  satisfy  the  symmetry 
It  is  shown  in  L 5 3 that  in  this  case  the  coefficients 
(36),  reduces  to 


4P  = V 

KL  Y 


df 


ds 


(43) 
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Hence,  the  normality  of  the  plastic  strain  rate  relative  to  loading  surfaces  in 
stress  space  holds  when  the  stress  response  is  specified  by  the  special  con- 
stitutive equation  (42). 


(D)  Materials  Whose  Stress  Response  is  Linear  in  (e^j-eP^) 
and  have  Constant  Coefficients 


Let  the  stress  response  be  specified  by  a special  case  of  (42)  in  which  the 
coefficients  LkLMTI  are  constants,  independent  of  e^Q  and  k.  Assuming  that  in 
the  initial  state  e^  = e^  = 0 when  sK^=0,  the  integration  of  (42)  yields 


<-  = l (p  -ey  ) 
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(44) 


where  the  constant  coefficients  LkXMN  satisfy  symmetry  conditions  similar  to 
(42)2* 

Inasmuch  as  the  stress  response  (44)  is  a special  case  of  (42)p,  it  follows 
that  the  normality  of  the  plastic  strain  rate  in  the  form  (43)  is  also  valid 
here.  However,  since  in  the  present  situation  the  material  coefficients  in 
(44)  are  independent  of  plastic  strain  and  the  work-hardening  parameter,  it  is 
possible  to  obtain  stronger  restrictions  than  those  given  by  (36)^  2 • Thus, 
by  starting  again  with  the  inequality  (25)  applied  to  a suitably  closed  spa- 
tially homogeneous  strain  cycle  (which  may  be  of  arbitrary  size),  we  obtain 
the  two  inequalities  [5] 


and 


(45) 
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where  s^j  in  (45)  refers  to  a state  of  stress  which  may  be  inside  or  on  the 
loading  surface  in  stress  6pace.  By  a procedure  similar  to  that  employed 
previously  by  Naghdi  [11],  from  the  inequality  (45)  we  can  prove  convexity  of 
the  loading  surface  f in  stress  space.  Hence,  as  long  as  the  coefficients 
^icLMN  in  (44)  are  constants,  the  work  assumption  (22)  for  the  strain  cycles 
also  implies  convexity  of  the  yield  surfaces  in  stress  space. 


(£')  Rigid- Plastic  Materials 

As  indicated  in  [5],  both  the  convexity  of  yield  surfaces  and  the  normality  of 
the  plastic  strain  rate  hold  in  this  case.  In  particular,  (36)^  can  be 
reduced  to 

V * ^ - v 5^  • ’s°  • <■"> 

which  is  a useful  expression  for  the  strain  rate  response  of  a rigid-plastic 
material.  By  suppressing  the  dependence  of  f on  eJJjj  and  K,  the  above  con- 
clusions hold  also  for  rigid-perfectly  plastic  materials. 

The  above  discussions  have  been  carried  out  so  far  in  the  context  of 
finite  deformation.  We  consider  now  one  further  case  concerning  the  small 
deformation  of  elastic-plastic  materials  whose  stress  response  is  of  the  form 
(44)  with  constant  coefficients. 


(F)  Elastic-Plastic  Materials  with  Small  Deformations 


To  avoid  the  introduction  of  additional  notations  we  continue  to  use  the  same 
notations  as  in  the  earlier  parts  of  the  paper,  even  though  in  the  linearized 
version  of  the  various  results  and  inequalities  the  distinction  between  the 
symmetric  Piola-Kirchhoff  stress  and  the  Cauchy  stress  disappears. 

W°  assume  that  the  stress  response  is  in  the  form  of  generalized  Hooke's 
law,  or  equivalently  of  the  form  (44)  with  the  infinitesimal  strain  e^  now 
taken  to  be  e^L=e^j  + ej^L,  ej^  being  the  infinitesimal  elastic  strain.  Then, 
the  conclusions  reached  previously  in  subsection  (D)  concerning  the  convexity 
of  the  yield  surfaces  and  the  normality  of  plastic  strain  rate  remain  valid 
here  also.  While  these  conclusions  are  the  same  as  those  which  follow  from 
Drucker's  postulate  [10],  the  consequence  of  the  linearized  versions  of  (29) 
and  (30)  are  not  the  same.  In  the  case  of  Drucker's  postulate,  we  have 


(ski.‘skl);£l  5 0 
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(48) 


While  the  combination  of  (43)  and  (45)  is  equivalent  to  (48)-^,  the  inequality 
(48)2  is  different  and  more  restrictive  than  the  linearized  version  of  (46)2* 
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